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We consider the diffusion scaling limit of the one-dimensional vicious walker model of Fisher 
and derive a system of nonintersecting Brownian motions. The spatial distribution of N particles is 
studied and it is described by use of the probability density function of eigenvalues of N x N Gaussian 
random matrices. The particle distribution depends on the ratio of the observation time t and the 
time interval T in which the nonintersecting condition is imposed. As t/T is going on from to 1, 
there occurs a transition of distribution, which is identified with the transition observed in the two- 
matrix model of Pandey and Mehta. Despite of the absence of matrix structure in the original vicious 
walker model, in the diffusion scaling limit, accumulation of contact repulsive interactions realizes 
the correlated distribution of eigenvalues in the multimatrix model as the particle distribution. 
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I. INTRODUCTION 



The vicious walker models, in which random walkers 
walk without intersecting with any others in a given time 
interval, were introduced by Michael Fisher and applica- 
tions of the models to various wetting and melting phe- 
nomena were described in his Boltzmann medal lecture 
Recently, using the standard one-to-one correspon- 
dence between walks and Young tableaux, Guttmann 
et and Krattenthaler et aZ.jjJ showed that exact 

formulae for total numbers of one-dimensional vicious 
walks, some of which were conjectured in previous papers 
jl], |^, ||, |(|, are derived following the theory of symmet- 
ric functions associated with Young diagram s ffy pj, ||] or 
the representation theory of classical groups |1C| ]. Impor- 
tant analogies between the ensembles of Young tableaux 
and those of Gaussian random matrices were reported by 
Johansson and then Baik |Q and Nagao and For- 
rester |L3|, |lj studied the vicious walker models using 
the random matrix theory |l5|, [l^] . 

The purpose of the present paper is to demonstrate 
more explicit relations among the vicious walker model, 
the symmetric function called the Schur function and the 
Gaussian ensembles of random matrices by considering 
the diffusion scaling limit of the one-dimensional vicious 
walks. Since each random walk converges to a Brownian 
motion in the scaling limit, the limit process of N vicious 
walkers will be a system of N nonintersecting Brownian 
motions fl7}| . In order to enumerate all possible nonin- 
tersecting paths of walkers realized on a spatio-temporal 



plane, we use the so-called Lindstrom-Gessel-Viennot for- 
mula 1 18, 191 E0|, which leads us to a useful determinantal 
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expression for the transition probability density of non- 
intersecting Brownian motions. We found that its initial- 
configuration dependence can be generally described by 
using the Schur function and the well-known properties 
of this function enable us to define the nonintersecting 
Brownian motions in which all particles start from a sin- 
gle position. Because the nonintersecting condition will 
be imposed for a given time interval, say T, all the par- 
ticles are immediately disunited from the initial point, 
and then they walk randomly keeping the nonintersect- 
ing condition. We have studied the time dependence of 
the spatial distribution of particle positions. We report 
in this paper that the position distribution of N nonin- 
tersecting Brownian motions can be identified with the 
distribution of eigenvalues of N x N complex hermitian 
matrix H coupled to a real symmetric matrix A, in which 
H and A are randomly chosen from the Gaussian ensem- 
bles. Such a two-matrix model was studied by Pandey 
and Mehta in which one parameter was intro- 

duced to control the coupling strength between two ma- 
trices. We will show that the time dependence of our 
process can be expressed by the parameter dependence 
of Pandey-Mehta's two-matrix model. 

Here we consider the probability density function of N 
real variables {x\, • • • , xn} with a real parameter /3 > 0, 

Pp(xi,---,x N ) = Ce-"XX/ 2 J] \ Xj -x k f 

l<j<k<N 

= Cexpi-pWdxj})) (1) 



with 
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j=l l<j<k<N 

where C is a normalization constant. It is known that (|l|) 
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with (3 = 1, 2 and 4 describe the distributions of eigenval- 
ues of random matrices in the Gaussian orthogonal, uni- 
tary and symplectic ensembles, respectively (abbreviated 
as GOE, GUE and GSE, respectively) fyf. For the one- 
dimensional TV Brownian motions, in which all particles 
start from the origin and the nonintersecting condition is 
imposed in the time interval (0, T], we will show that (i) 
at the very early stage, i.e. i/T < 1, the particle distri- 
bution is described by using GUE, (ii) as time t is going 
on, a transition from GUE to GOE is observed, and (iii) 
at the final stage t = T the particle distribution can be 
identified with GOE. As shown by the second equality 
of (Q), the Gaussian ensemble of random matrices can 
be regarded as the thermodynamical equilibrium of one 
dimensional gas system with (two-dimensional) Coulomb 
repulsive potential (||) at the inverse temperature [3. Here 
it should be noted that the vicious walkers on a lattice 
have only contact repulsive interactions to satisfy the 
nonintersecting condition. The global effective interac- 
tions among walkers are accumulated by taking the diffu- 
sion scaling limit and as its result a long-ranged Coulomb 
gas system is constructed. Such emergence of long-range 
effects in macroscopic scales from systems having only 
short-ranged microscopic interactions is found only at 
critical points in thermodynamical equilibrium systems, 
but it is a typical phenomenon observed in a various in- 
teracting particle systems in far from equilibrium. 

In particular, in the limit T — > oo, that is, when the 
nonintersecting condition will be imposed forever, we can 
derive a system of stochastic differential equations for the 
process with the drift terms which act as the repulsive 
two-body forces proportional to the inverse of distances 
between particles. In other words, the scaling limit of 
vicious walks with T — > oo can realize Dyson's Brownian 
motion model at [3 = 2 p3| . It is reasonable to obtain 
such a stochastic process from the vicious walker model, 
since it is known that Dyson's Brownian motion model 
at f3 = 2 can be mapped to the free fermion model JlT], 
Izjf . The transition from GUE to GOE is, however, first 
reported for vicious walkers with T < oo and explained 
using the two-matrix model in the present paper. 



II. MODEL AND 
LINDSTROM-GESSEL-VIENNOT 
DETERMINANT 

One-dimensional vicious walks are defined as a subset 
of simple random walks as follows. Let {i?^. J }fc>o ; J — 
1, 2, • • ■ , N, be the N independent symmetric simple ran- 
dom walks on Z = {• • ■ , —2, — 1, 0, 1, 2 • ■ •} started from 
N distinct positions, 2s\ < 2s 2 < ■•• < 2sn, sj G Z. 
That is, 



Rn 



2Sj 



and R 



k+1 



Rl 



1 or R' + 1, 



for j = 1, 2, • • • , N, k = 0, 1, 2, • • • . Fix the time interval 
K as a positive even number. The total number of walks 
is 2 NK , all of which are assumed to be realized with equal 



probability 2~ NK . We consider a subset of walks such 
that any of walkers does not meet other walkers up to 
time K . In other words, the condition 



Rl 1 <R{ 2 < ■ ■■ < Rl N , 



k = 1,2, 



(3) 



is imposed. Such a subset of walks is called the vicious 
walks (up to time K) @, g. Let N N (K; {ej}\{sj}) be 
the total number of the vicious walks, in which the N 
walkers starting from 2s\ < 2s 2 < ••• < 2sn arrive 
at the positions 2ei < 2e 2 < ••• < 2e^ at time K. 
Then the probability that such vicious walks with those 
fixed end-points are realized in all possible random walks 
started from the given initial positions, which is denoted 
as V JV ({i#}£ =0 ;J22 = 2e J ), is 



V N{{Rk}k=0\ R K = 2e J 
We also consider the probability 



N N {K-{e 3 }\{ Sj }) 



yNK 



VnUR^Lo) 



E 

ei<e2<---<ejv 



MK'}£=o;*2 = 2e 3 -). 



Consider a subset of the square lattice Z 2 , 

Ck = {(x, y) G Z 2 : x + y = even, < y < K}, 

and the set £k of all edges which connect the nearest- 
neighbor pairs of vertices in Ck- The lattice (Ck,£k) 
provides the spatio-temporal plane and each walk of 
the j-th walker, j = 1, 2, • • • , N, can be represented 
as a sequence of successive edges connecting vertices 
Sj = (2sj,0) and Ej = (2ej,K) on it, which we call 
the lattice path running from Sj to Ej. If such lattice 
paths share a common vertex, they are said to inter- 
sect. Under the vicious walk condition (|J), what we 
consider is a set of all TV-tuples of nonintersecting paths 
fj"of . Let 7r(S — > E) be the set of all lattice paths from 
S to E, and 7r ({S J -}jl 1 -> {Ej}^) be the set of all 
iV-tuples (iti, ■ ■ ■ , 7r/v) of nonintersecting lattice paths, in 
which 7Tj runs from Sj to Ej, j = 1,2, •••,7V. If we 
write the number of elements in a set A as \A\, then 

N N (K;{ ej }\{^}) = ko({S,}f = i - {Ej}f =1 )|. 

The Lindstrom-Gessel-Viennot theorem gives [[l8| [l|, 
H (see also i 10), 

N N (K;{ ej }\{ Sj }) = det (|7r(S* -> Ej)|). 

1<J, k<N 

Since |7r(Sfe — > Ej)| = (#-/2+aj,-e-)' we have the following 
binomial dctcrminantal expressions 



VN({R s k i }^ ;R s i = , 2e j ) 



-NK 



det 



K 



i<j,k<N \ \K/2 + s k — ej 



, (4) 



and 



2 -NK V" 



ei<e2<-<ejv 



det 

l<j,k<N 



K 

K/2 + s k 



(5) 
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III. SCALING LIMIT OF VICISOUS WALKS 

Recently Krattenthaler et al. || evaluated the asymp- 
totes of (0) for large K in the two special initial- 
configurations, (i) Sj = j — 1 and (ii) Sj = 2(j — f), 
as 

W({i$'}f=o) = a N b N ({ Sj })K- N ^ N -^ (1 + 0{l/K)) , 

(6) 

where 

N/2 

(2 N /ir) N / i Y[(2j - 2)! if N = even 

"■^ ^ J_1 (iV-i)/2 ( 7 ) 

(2 w+i /7r) (JV-i)/4 TJ ( 2 j-l)! if TV = odd, 

3=1 

and 

Mti - 1}) = 1. b N ({2(j 1)}) = 2 Ar ( JV " 1 )/ 2 . (8) 

We found that their result can be immediately gener- 
alized as 

b N ({s(j 1)}) = s N ( N -^ 

for Sj = s(j— 1), s = 1,2,3,- • • . This observation suggests 
that we can take the scaling limit L — > oo, where the time 
interval K oc L and the initial spacing of walkers s oc sfL. 



Therefore it is concluded that the ratio of two determi- 
nants A\/ Aq is a polynomial in zx, ■ ■ ■ , z/y. Moreover, 
it can be readily seen from (||) that the Schur func- 
tion is a homogeneous polynomial of degree JZj—i Xj m 

Zl,---,Z N . 

Let q be a complex variable and set Zj = g* 3 ' - 1 in (|^). 
Then we have 



sa(1,<Z,<7 2 , 



•••,9 



det ( g (j'-i)(>*+^-*)) 
det (qU-mN-k^ ■ 

l<j,k<N 



Appropriate application of the formula of Vandermonde 
determinant ( ^p| ) gives the product form 



q 



jXj—Xk+k—j _ ^ 
9 fe "J - 1 ' 



l<j<fc<AT 

Taking the limit q — > 1, we have the formula 

Xj - Xk + k- j 



SA (i,v--,i)= n : ^ 

l<j<k<N 



k-j 



(11) 



The Schur function is a character of the irreducible rep- 
resentation specified by A of the group GL(V) and Jll] ) 
gives the dimension of the representation. 



B. Diffusion scaling limit 



A. Schur function 

In order to describe the scaling limit of the vicious 
walks, the symmetric function called the Schur function is 
useful. Here we give some of the fundamental properties 
of Schur function |?], S, |9|, which will be used below. 

A partition A = (Ai, A2, • • • , Ajv) is a non-increasing 
series of non-negative integers, Ai>A 2 >--->Aat>0. 
Let V be the ./V-dimensional complex vector space. Then 
the Schur function s\(zi, ■ • • , Zn) associated with A is a 
function of (21, • • • , zn) G V defined by 



s\(zi, ■ ■ ■ ,z N ) 



det 

l<j,k<N 



( X k +N-k\ 



det (zf- k ) 

l<j,k<N y J ' 



(9) 



Let A\({zj}) be the numerator of (Q), which is an TV x N 
determinant. If we set zi 1 — zt 2 for 1 < £\ < £2 < N, 
then A\({zj}) = 0, since the £i-th row is equal to the 
^2-th row. Then it is divisible by each of the differences 
Hi — zg 2 ,l < £1 < £2 < N, and hence by their product 
rii<j<fe<Af( 2: i — z k)- This product of all differences is 
known as the Vandermonde determinant, which is noth- 
ing but the denominator of (||); 

A ({zj}) = del : (*?-")= II (^-^). (10) 

l<1,k<N J 

- l<j<k<N 



We set 



K = Lt, Sj = 



VI 



vT 



"3 ' 



(12) 



for j = 1, 2, ■ • • , N, and take the limit L — > 00. Since in 
this limit each random walk R s k 3 converges to a Brownian 
motion, whose distribution function solves the diffusion 
equation, this scaling limit is especially called diffusion 
scaling limit. First we remark that, for each strictly in- 
creasing series of integers yi < 1/2 < ■ ■ ■ < un, a weakly 
decreasing series of integers £(y) = (£1(2/), ■ • • , £n(v)) can 
be assigned by setting 



£j(y) = y N -j+i - (JV - j), j = 1,2, • • • ,N. 



(13) 



Then we can prove that, for given t > 0, x\ < X2 < ■ ■ ■ < 
x N , and yi < y 2 < ■ ■ ■ < vn, 



(r~\ N 



LI 



k=0 



■R^=VZ V j 



(2wt)- N / 2 det 

l<j,k<N 



exp 



2t 



(Xk -Vj)' 



-N/2, 



,X\/t x 2 /t _ _ _ x N /t 



x exp 



(-hp 



x*+y2) h N {{e'>l*}), 



(14) 
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where s^t y -\(zi, • • • , Zn) is the Schur function associated 
with £(y), defined by (||) with A = £,(y), and 

h N {{z 3 }) = det (z*- 1 ) 

l<j,k<N J 

= (-l) N ^ 2 A ({ Zj }) 

= n i^-zj). (15) 

l<j<k<N 



The proof is given as follows. Setting (12), we apply Stir- 
ling's formula to the RHS of (§) multiplied by (\/I/2) JV , 



lim 2-^* I ^ 


L — >oc 


V 


det 


lim 2" 


l<j,k<N 


l L — ► oo 


det 


l 

■ e 


l<j,fe<JV 


.%/27rf 



det 



Lt 



<N \ \ — 4- ^( X k-Vj) 
V 2 2 



"(f) 

■(x k -yj) 2 /2t 



Lt 



(16) 



which gives the first equality of (fbil). For the second 
equality, we rewrite ( |T6| ) as 

(2f()-" /2 e^ ( ^'' )/2i det f e W*V 

l<j,/c<Ar V / 

The determinant is written as 
det f e*^/*) 

l<j,k<N \ J 

det 1 < J , fc < jV ((e^/*)^-^ 1 ) / 
~ /7 x /t\N-k\ h N{{e })■ 

Using (|l^) and the definition of Schur function (||), the 
second equality of (|l4| ) is obtained. 

We consider the rescaled one-dimensional lattice 
Z/(\/Z/2), where the unit length is 2/yL, and let R% 
denote the symmetric simple random walk starting from 
x on Z/(vT/2). Then (fj) implies that 

km U W ({$% }f n ; ^ € [%, % + d% 

L— >oo y L J fe=0 

= fN&{Vi\\{Xj\)d N V . 

Here we can give two expressions for /iv(i; 

/jv(*;{yi}|{^i}) 

= (2^t)- Ar / 2 det 

l<j,k<N 



exp ( -^(^fc - 



f-^E^+^^c^/*}). a?) 



x exp 



Since the vicious walkers are defined by imposing the 
nonintersecting condition (||) up to a given time K, the 



process depends on the choice of K. That is, the process 
is temporally inhomogeneous. This feature should be in- 
herited in the process obtained in the diffusion scaling 
limit. Since each random walk converges to a Brownian 
motion in the diffusion scaling limit, the limit process of 
the TV vicious walkers can be called the N nonintersect- 
ing Brownian motions. Let T > and we consider the 
iV nonintersecting Brownian motions in the time interval 
(Q,T]. Set 



Af N (t;{ Xj }) 



d N y f N (t;{ yj }\{ Xj }). 



Vi<—<yN 



For < s < t < T, x\ <•••< xn,Vi < ■ ■ ■ < yN, 
the transition probability density from the configuration 
{xj} at time s to {yj} at t is given by 



a T (s ix Vt iv \) - Mt~s;{ yj }\{^}WN(T-t;{y 3 }) 

(18) 

since the numerator in RHS gives the nonintersecting 
probability for (0, T] specified with the configurations 
{xj} and {yj} at times s and t, respectively, and the de- 
nominator gives the probability only specified with { X j} 
at s, where we have used the Markov property of the pro- 
cess. The temporal inhomogeneity is obvious, since RHS 
depends not only t — s but also T — s and T — t. 



C. t — » oo asymptote of A/jv(i; { x j}) 

It should be noted that, since A/jv(t; {xj}) is the in- 
tegral of /jv(t; over all possible end-positions 
{yj}, it is the probability that N Brownian motions start- 
ing from {xj} do not intersect up to time t. Before study- 
ing the stochastic process defined by the transition prob- 
ability density (|l8|), here we assume |x| = YljLi \ x j\ < 00 
and evaluate the t — > oo asymptote of A/jy(*; { x j})- I n or " 
der to do that, the second expression of /jv(i; {yjlK^}) 
in (Rt|) will be useful, 



■h N ({e x ^}) 



{2nt) N / 2 

d N y s i(y) {e^\--- 1 e^/ t )e-Y.yV 2 K 



iy\<—<VN 
By @, © and ©, 



lim^ (j/) (e-/*, 



)=*ec*)(M, 

m{%»/ n (fc-i), 

l<j<k<N 



and 



t — >OC 

We define 



li m t Jv ^- 1 )/ 2 / lw ({e^ /t }) = %({^}). 



= ^({^})/ n ( fc -j)- 

l<j<k<N 
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Remark that this definition of 6at({x}) is consistent with 
(H). Then we have 



M N {t;{x 3 })=t- 



N 2 /2 t>N({Xj}) 



(2ir) N / 2 Nl 

x / d N y e-T,vV»\h N ({ yj })\ x (1 + 0(1/*)) 



f -N(N-l)/4 fyv({Zj}) 
( 27r )JV/27V! 



d^M e-E^/ 2 |/ lAr ({ % })| X (1 + 0(1/*)) 



as * tends to infinity, where we have used the facts 
that with the absolute values the product of differences 
|/iAr({t/j})| is invariant under permutation of yj, and 



yj/Vi. The last integral is the special case (7 = 1/2 



and a = 1/2) of 

Jd N ue- a ^ J] 



l<j<k<N 



N 



(2ir) N / 2 (2a)- N ^ N -V+V/ 2 TT ^ ( / 1+j7 ^ (19) 

J -- L r(i + 7) 



i=i 



which is found in Mehta ||l5f ( eq.(17.6.7) on page 354), 
whose proof was given in |25| by use of Selberg's in- 
tegral (26). Here T(x) is the Gamma function with 
the values T(3/2) = ^/2 and l\f=i r (! + J'/ 2 ) = 

2 -N(N-l)/2^i^N m fljVj where aN ig giyen by ^) 

Then we have 

Af N (t;{xj}) 

= t^' N 2- 2 ^a N b N ({x 3 })(l + 0(l/t)) (20) 



with 



1> 



N 



-N(N-l), 



(21) 



as t tends to infinity, where ^at is known as the critical 
exponent of survival probability of vicious walkers jl], ^, 
|27|, ||. Since 

t-^2-W»b N ({ Xj }) = (Lt)-^b N ({VL Xj /2}), 

© suggests that the result (|) with (@) and (|) of Krat- 
tenthaler et al. shall be generalized for arbitrary initial 
positions of vicious walkers on the lattice. 



IV. GAUSSIAN RANDOM MATRIX 
ENSEMBLES AND DYSON'S BROWNIAN 
MOTIONS 

In this section we study two special choices of T ; T = t 
and T — ► 00. We show that there is an interesting corre- 
spondence between these choices of T and the Gaussian 
ensembles of random matrices. In order to see it we con- 
sider the limit |x| — > 0, where |x| = Ylj=i \ x j\- ^ wm be 
shown that the second expression of fN(t]{yj}\{%j}) hi 
( |l7| ) is useful for taking this limit. 



A. T = t case and GOE 



Since the first expression in ( |l7| ) gives 

lim t ^o/jv (*;{%•} I {xj}) = I\f=i s ( x j - Vj) with 
Dirac's delta functions, A/jy (0; jz 7 |) = 1 for any {xj}. 
Then setting T = t makes (|l^) depend only on t — s. 
Set s — and use the second expression in (|l7]) for 
/jv(£; {yj}\jxj}) and A/at(*; {xj}). By virtue of the Schur 
function (pd|), for < > and |x| <C 1, we have 

x S5(y) (i,...,i) e -E^ 2 / 2 *(i + o( N )) 



(271)^/2 



^ 2/2t M{%}) 



n 



l<j<k<N 



k-j 



x (i + o(H)), 



and 



Af N (t;{xj}) = {2Tit)- N ' 2 h N ({e x ^ t }) 

d N y s m (l,---,l)e-^ /2t (l + 0(\x\)) 

e Xk/t _ e Xj/t 



n 



yi<---<VN 
±N{N-l)/i 



(2tv) n / 2 c n 11 k-j 

y 1 " l<j<k<N J 



where the integral ( |19| ) was used and 

2N(N-2)/2 



(i + o(N)), 



TV 



Cat 



2 W/2 II r 0V2) 



Then (||) gives 

^(0,{0};t,{ yi }) = cjvt^e-E^ 2 / 2 *^^.}) 
for yi < ■ ■ ■ < vn with 

Cat = ±N(N + 1). 

It means that 

g t N (0,{0};t,{y j }) = N) <?£ OB (M;*) 
for yi < ■ ■ ■ < ijn , where 

^° E ({2/,};- 2 ) = ^ 2CN exp f-^f»j Mfoi}) 

(22) 

is the probability density function of eigenvalues {yj} of 
random matrices in the Gaussian orthogonal ensemble 
with variance a 2 ITHl- 
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T -> oo limit and GUE 



C. Dyson's Brownian motions 



Let 



p N (s,{xj};t,{yj}) = Jim g£(s, {:£,-};*, {%}) 



By use of Q2^) we can determine the explicit form for any 
initial configuration x± < X2 < ■ ■ ■ < xm in this case as 



p N (0,{xj};t,{yj}) 



h N ({xj}) 



Mt;{ yj }\{*j})> ( 23 ) 



where Hm is given by (|15|). Moreover, if we take the limit 
|x| — * 0, we have 

p N (P,{0YMVi}) = c' N t-^e-^-/ 2t h N ({ yj }f, (24) 
with 



N 2 

Cjv = — and c ' 



JV 



JV 



That is, we have the identity 

p N (0,{0};t,{ yj }) = N\ g% UE ({ yi };t), 
for f/i < • • • < t/jv, where 



<& uii (te};0 



-^■er 2Cn exp 
AH 



(25) 

is the probability density function of eigenvalues {j/?} of 
random matrices in the Gaussian unitary ensemble with 
variance cr 2 p5[ . 

In the case T — > oo, the nonintersecting condition will 
be imposed forever, while in the case T = t, there will be 
no condition in the future. The distributions of particles 
at present depend on the condition in the future. 

By generalizing the calculation, which we did in the 
case T = t, for arbitrary T and comparing the result 
with (p4h, we have 



g%(0,{0};t,{Vj}) = - T *„ MN(T-t;{ yj }) 

P N(p,{oYMvi}) ° N h N ({ yj }) 



(26) 



for yi < ■ ■ ■ < j/jv, with (Elf) and 



Cjv 



CJV 



'JV 



JV 



x N / 2 TT 

/i r(i/2) 



When jV = 2, we can consider the process of one variable 
y = 2/2 — 2/1 • I n this case g% and ^2 define the Brow- 
nian meander and the Bessel process, respectively, both 
of which are stochastic processes well-studied in proba- 
bility theory ^9|. The equality (^6|) can be regarded as 
the multi-variable generalization of Imhof 's relation |3(J 
between the Brownian meander and the Bessel process. 



In the limit T — ► oo we have obtained the compact ex- 
pression ( pi) ) for any x\ < ■ ■ ■ < x^ and J/i < • • • < J/jv- 
In this section, we show that a system of stochastic differ- 
ential equations can be explicitly derived for fl23| ) . Using 
it we will explain why we have the GUE distribution. 

Let 



JV 



Ek({x j })= Yl 
It is easy to verify that 



1 



for k = 1,2,- --N. 



E k{{Xj}) = —logflNdXj}), 



for k = 1,2,- ••, N, and 

JV r 



E 

k=l 



A- Ek{{Xj}) + {Ek{{Xj}) f 



(27) 



0. (28) 



Using these equalities, we can prove that 
Pjv(0, {xj}] t, {yj}) solves the equation 

q I N d 

-u{t-{x 3 }) = -AufaixjV+^Ekiix^—ufcixj}) 



fe=i 



dx k 



(29) 

where A = J2k=i ® 2 ' l® x \- The proof is the following. 
First we remark that the first expression in (|lj]) states 
that /at is a finite summation of the products of Gaussian 
kernels and thus it satisfies the diffusion equation 
Therefore 

|p.(0,{* J } ; t > {y J }) = i^g^(t;{y i }|{ a!i }). 

Then we can find that, if {Ek{{xj})} satisfy the equa- 
tions 



N 1 







fe=l 



JV 

E 



9 1 1 f 9 



/jv(t;{j/j}|{a:j}) 

M*;{i/i}|{*i})K 3 0) 



and 



(9 



fe=i 



5a; fe /iAr({xj}) 



1 w 

IE 



<9 2 



2 dx 2 k h N ({xj})'- 



(31) 

( P9| ) holds with u (t; { xj}) — pjy(0, {xA; t, {y 7 |). It is 
easy to see that ( |30| ) is satisfied if (|27[) holds for any 
k = 1,2,---,jV. Moreover, using (p7[), we can reduce 
( pl| ) to (|28|). Then the proof is completed. 

The above result implies that the process defined in the 
limit T — > oo is the system of N particles with positions 
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x%(t), X2(t), ■ ■ ■ ,xpf(t) at time t on the real axis, whose 
time evolution is governed by the stochastic differential 
equations 



dx k {t) = E k {{ Xj (t)})dt + dB k (t), 



(32) 



k = 1,2, ■ ■ ■ , N, where {Bk(t)}^ =1 are the independent 
standard Brownian motions; 

B i (0) = 0, (B j (t)) = 0, 

((B 3 (t) - B,{ S )){B k {t) - B k (s))} =\t- s\S jk 

for any t, s > 0, j, k = 1,2, ■ ■ ■ , N . Because of the scaling 
property of Brownian motion, \faBj (t) is equal to Bj (at) 
in distribution for arbitrary a > 0. Then, if we set t = 2t' 
and write x k {t) = Xk(t'), (|32|) are in the a = 0, (i = 2 
case of the equations 





dx k {t') 
k = l,2,---,N, with 



Pa=-W a ({xj(ir)})<tf + V2 dBfc(i'), (33) 
ox k 



n ~2 



l<j<k<n 



When a = 1, (|33| ) is known as the stochastic differ- 
ential equations for the Dyson Brownian motions at 
the inverse temperature (3 and the stationary distribu- 



tion ex exp(— /3H /1 ({i J })) 



If a = 0, the factor 



cxp(— /3a will be replaced by exp(— Xj/4±') 

for finite £' and thus when t' = a 2 / (3 we may have the 
Gaussian distribution oc exp(— (/3/4cr 2 ) £|)ftjv({ij})^- 
Setting (3 = 2 gives the form (|25|). 

It should be noted that the system of diffusion equa- 
tions describing the Dyson Brownian motions with j3 = 2 
can be mapped to the free fermion model jl7], Q . 

For general T < oo, we will have the stochastic differ- 
ential equations 



be identified with the eigenvalue distribution of random 
matrices in GOE, and (ii) if T — * oo, it becomes GUE. 

For a fixed T < oo, the above results are summarized 
as follows. Consider the one dimensional N Brownian 
motions all starting from the origin at time t = 0. We 
impose the nonintersecting condition for the time interval 
(0,T]. As the ratio t/T —> 0, the particle distribution is 
asymptotically described by GUE. On the other hand, at 
t = T, it can be identified with GOE. This implies that as 
time t is going on from to T, there occurs a transition 
of distribution from GUE to GOE. In this section, we 
study this transition. 

GUE is the ensemble of complex hermitian matrices 
and GOE is that of real symmetric matrices. The de- 
grees of freedom are, when the matrix sizes are N, N 2 
and N(N + 1)/2, respectively If we change the variables 
from these independent matrix elements to the eigenval- 
ues and other mutually independent variables, and then 
if we integrate the distribution functions over all vari- 
ables other than eigenvalues, we will have the probabi lity 
density functions for TV real eigenvalues as (g^) and ( [22] ) 

Although the vicious walker model has no matrix 
structure at all, here we show that its diffusion scaling 
limit, nonintersecting Brownian motions, can be regarded 
as the reduction of a one-parameter family of ensembles 
of matrix structures to a variable space of eigenvalues. 
The "hidden structure" is not a single matrix but a two- 
matrix model, in which a complex hermitian matrix is 
coupled with a real symmetric matrix. 

In the first subsection we will derive the two-matrix 
model from the nonintersecting Brownian motions and 
the transition from GUE to GOE will be discussed in 
the second subsection. In the third subsection we will 
show that the obtained two-matrix model can be iden- 
tified with the two-matrix model of Pandey and Mehta 
PH ^2) by appropriate scale transformation of matrix 
elements. 



dx k {t) = El{{x {t)})dt + dB k (t), 



A. From vicious walker model to two- matrix model 



for k = 1,2,- • •, N, with 



El{{xj}) 



d 

dxk 



lnAf N {T-t;{ Xj }). 



TWO-MATRIX MODEL 



In Section III we have constructed a system of non- 



intersecting Brownian motions in one dimension as the 
diffusion scaling limit of vicious walks. The obtained 
transition probability density (|l^) is temporally inhomo- 
geneous and the particle distribution depends not only 
the observation time t — s but also on the time inter- 
val T, in which nonintersecting condition is imposed. In 
the case that all particles start from the origin at time 
s = 0, it was shown in Section [V that, (i) if T = t, it can 



The generalized Imhof relation (|26| ) with ( f24|) gives 
g%(0, {0}; t,{ yj }) 
oc e-T.vV 2 % N ({y 3 }) J d N zsgn(h N ({ Zj })) 

1 



x dct 

l<j,k<N 



exp 



2(r_t)^ 3 

where sgn(x) = The RHS is rewritten as 

ft JV ({%■}) / d N z sgn(h N ({zj})) 



(34) 



x dct 

l<j,k<N 



1 

2? 



(Vj ~ z k)' 



2(T-t) 

= h N {{y 3 }) J d N z S gn(h N ({ Zj }))e-J: z ^ 2T 
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x det 

l<j,k<N 



exp ( I . 



Setting (tjT)zj — dj, j = 1, 2, • • • , N, we have 
<?£(0,{0};t,{ % }) 



JV 



X /'.\ ({'/,!-') / {<!.,•})) exp I T^lX/'./ 

T 



x det 

l<j,k<N 



exp 



2t(T - i) 



(35) 



Consider an ensemble of N x AT real symmetric matrices 
{A} with an integration measure 

dA = J] dA jk . 

l<j<k<N 

Let {ai, a2, • • • , a/v} be the eigenvalues of the matrix A 
and {pi,j?2, ■ • • ,Pjv(/y-i)/2} be other mutually indepen- 
dent variables. Then 

N N(N-l)/2 

dA = J({aj}, {pj}) Y[ da 3 Y[ d Pki 
3=1 k=l 

where J({aj}, {pj}) is the Jacobian 

d(An,A 12 , ■ ■ ■ , A NN ) 



J{{°>j},{Pj}) = 



d(au ■ ■ ■ , on, Pi, ■ ■ ■ ,Pn(n-i)/s) 

It is known that we can write 

J({a j },{Pk}) = \h N ({a j })\f({ Pk }), 

where f({pk}) is independent of dj's (jl5). Therefore, 
for any function G({cij}) of {a\, ■ ■ ■ , a^}, we have the 
identity 

N 

dA G({a 3 }) =c/n doj l^({%-})|G({a,-}) (36) 
J 3=1 



with 
Set 



, N(N-l)/2 



[ n /({Pi»- 
fc=i 



TV 



I IT 

sgn(/i JV ({a J }))exp -— 



IMW)! 



x det 

l<j,k<N 



2t 2 ^ 3 

3 = 1 



exp 



r 



2t(T-t) 



(yj ~a k )'- 



Then using the formula (p6|), (35) becomes 
gjr(0, {0};t, { Vj }) 



oc h N ({yj}) / gL4 



/l/v({Oj}) 



exp 



x det 

l<j±<N 



exp 



2t(T - t) 




(Vj - a k) 



Next we use the following integral formula |3^, [33|, [34j ; 
for x N hermitian matrices A and B having eigenvalues 
{ai, • • • , ajv} and {61, • • • , &a?}, respectively, and for any 
constant 7, 



J dU exp [7 tr(A - U^BU) 2 ] 



1 



h N({aj})h N ({bj}) i<j,k<N 

where the integral is taken over the group of unitary ma- 
trices U. Then ( |37| ) can be written as 

ff^(0,{0};t,{y i }) 

« Mfe}) 2 j dU j dA exp f-^trA 2 



x exp 



2t(T - 4) 



tr(C/ t y"C7 - A) 2 ^j , 



(38) 



where Y is the N x N diagonal matrix such that Yj k — 
VjSjk- Since U is a unitary matrix, H = U^YU is an 
N x N complex hermitian matrix. Then the integrand 
of ([38]) can be regarded as a weight for two matrices H 
and A given as 



exp (-tr (j H H 2 - jhaHA + j A A 2 )) 



with 

1H = 



T 



2t{T - t) 



, 7ha = 



t(T - t) 



, 1A 



2t 2 (T-t)' 



(39) 

Consider an ensemble oi N x N complex hermitian 
matrices {H} with the integration measure 

dH= Y[ dRe(H jk ) 11 dlm(H jk ). 

l<j<k<N l<j<k<N 

For each complex hermitian matrix H, let {j/i, • • • , yjv} 
be a set of eigenvalues and U be the NxN unitary matrix 
such that H = U'YU with Yjk = VjSjk- Then it is known 
that the integration measure dH can be factorized into 
the product of the Haar measure for unitary matrices dU 
and an integration measure for eigenvalues [nj 

JV 

dH cx dU x h N ({y 3 }) 2 dy 3 . 

3=1 

Now we introduce a two-matrix model, which consists 
oi an N x N real symmetric matrix A and an N x N 
complex hermitian matrix H, with a probability density 
function 

fi N (H, A) = — exp (-tr( lH H 2 - j HA HA + j A A 2 )) . 



Here ^h^IhAiIa are given as ( |39j ) and Zn is the parti- 
tion function of the two-matrix model, 

Z N = j dH j dA exp (-tv{ lH H 2 - lHA HA + lA A 2 )) . 
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Then the relation 

9 %(0,{0};t,{ yj }) <x h N ({ yj }) 2 J dU J dA m (U^YU,A) 
is established. 

B. Transition from GUE to GOE 

Consider the Gaussian ensembles of real symmetric 
matrices {A} and complex hermitian matrices {H} with 
sizes N with the probability density functions 

v N {A) = CUexp (_— T trA s 

and 



v N (H) = C H exp 



2a% 



tvH' 



respectively, where 

and C A = 2-^/2(^2 )-<»,C H = 2"^ (7^)"^ . Then 
consider the convolution 



fa(H) = J dA v N {A)v N {H - A). 

Since, for 1 < j, k < N, 

H jk =Re(H jk ) + iIm{H jk ) 
with i = \J — 1, and 

Re(A jk )=A jk , Tm(A jk ) = 0, 
the convolution is also Gaussian distribution in the form 

( (Re(H jk )) 2 , (Im^-fe)) 5 



oc exp 



\ -< I Kf^Mjkl 



I) 2< 



Then (pf) gives 



gir(P,{OYMVj}) « Mfe}) 2 / ^ Aiv(-ff) 



(Re^-fe)) 2 , (Im(^- fc )) 2 



2a T 2 

In 



where H = WYU and 



(40) 



(41) 



Now the transition from GUE to GOE is explicitly 
represented by the time-dependent variances ©. With 
a fixed finite T, if < t < T, cx^ = i ~ cr 2 m . Then the 
real and imaginary parts of complex hermitian matrix 
elements are equally distributed as in GUE. While a^ c 
increases linearly in t, <7 2 m increases in time t only up 
to time t = T/2 and then decreases in time. At time 
t = T ', a -2 = 0, which implies that the imaginary parts 
of matrix elements are zeros with probability one. Then 
the distribution is identified with GOE. 



C. Pandey-Mehta's two- matrix model 

As an interpolation between GUE and GOE, Pandey 
and Mehta introduced a family of Gaussian ensembles of 
hermitian matrices {H} with one parameter a G [0, 1] 




(MH lk ) f , (Im(g jfc )) 2 ' 



Av 2 



where 



= {2(1 



4v 2 a 2 j 
and Cpm 



,(42) 



2 -N/2 a -N(N-l)/2 ( 2 7TV ^-N*/^ get 



<(2T - t) 



T 



Then, it is easy to see that, if 



the equality 



a 2 = 1 , 

T 



K "fi N ( K H) = ^ N M (H,a) 



(43) 



(44) 



(45) 



is established. 

For an even integer N and an antisymmetric N x N 
matrix B = (bj k ) we put 

Ph<j<k<N{bjk) 
= (N/2)l ^ S S n ( (T )^(l)^(2) & <r(3) -(4) ■ ■ • &<t(JV-1>(JV), 

where the summation is extended over all permutations 
cr of (1,2,..., N) with restriction a(2k - 1) < a(2A;), 
fc = 1, 2 , .... N/2. This expression is known as the Pfaf- 
fian |M . Pandey and Mehta showed that the probability 
density function of eigenvalues {j/j} of the complex her- 
mitian matrices, which are distributed following (0), is 
given by 



9n = Cat (a) exp 



x/ ) , Ar ({y J })Pf 1 < J<fe < A r(F :;fe ), (46) 



10 



where, setting 



l-a 4 



1/2 



exp 



l-a 4 
4a 2 



V 2 dy, 



if N is even, F jk = f(y j - Vk),j,k = 1,2,- • -,N, 
and if N is odd, we use above and in addition 
Fj,N+i = —Fn+ij = l>i = 1) 2, • • • , N, Fn+i,n+i — 
0, and CN(a)- 1 = 2 3JV / 2 (1 - a 2 ) N( - N ^/ 4 (l + 
a ^)-N(N+i)/4Yl^ =i T(l + j/2) H §§. Then the rela- 
tion (^) and the equality (|||) with (|43|) and (j||) imply 
the expression 

<?£(0, {0}; t, {%}) = c^T JV ( w - 1 )/ 4 t- w2 /2 ftjv({2/ . }) 



x exp 



(47) 



where 



iV 



AT 



if A" is even 
if TV is odd, 



and 



Fj k {t,{yt}) 



=Erf 



2/fc ~ % 
2V< 



1 

-1 




if 1 < j, k < N, 



if 



1 < j < N,k = N + l, 
if j = N + 1,1 < k < N, 
if j — k — N + 1 



with 



which the nonintersecting condition is imposed, and it 
can be described by use of the probability density func- 
tion of eigenvalues of N x N random matrices in Gaussian 
ensembles. It was shown that the particle distribution 
depends on the ratio t/T and a transition from GUE dis- 
tribution to GOE distribution occurs in its time devel- 
opment. Such a transition between different ensembles 
of random matrices had been studied in the two-matrix 
model by Pandey and Mehta, in which a parameter a was 
introduced and a one-parameter family of random ma- 
trix ensembles was considered. The present work showed 
that the scaling limit of vicious walk model realizes such 
a two-matrix model as a stochastic process, in such a way 
that the parameter a is continuously changed following 
(0) as the system is developing in time t up to T. 

In the present paper, since we have considered only 
the transition probability density between two different 
times, a two-matrix model was analyzed. As briefly re- 
ported in [ j36| , however, multi-time correlations among 
particles at intermediate times between t = and 
t = T can be identified with inter-matrix correlations 
in a multimatrix model. The corresponding multima- 
trix model for calculating M-intermediate-time correla- 
tion functions is a version of Nagao's (M + 1) matrix 
model, where one real symmetric matrix is combined at 
the end of a chain of M complex hermitian matrices j[7| . 
This observation implies that the diffusion scaling limit 
of the vicious walker model is mathematically identified 
with a matrix chain, which is set along the time axis. 
The time development from GUE to GOE in the process 
can be then regarded as appearance of an edge effect, as 
the observation time t on the time axis is approaching 
to the end point t = T in this chain structure. Further 
study on the relations between multimatrix models and 
nonequilibrium interacting particle systems is desired. 



Erf(x) 



du e 
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APPENDIX: INTEGRATION VERSION OF 
OKADA'S MINOR-SUMMATION FORMULA 



VI. CONCLUDING REMARKS 

In the present paper we performed the diffusion scaling 
limit of vicious walker model in one dimension and con- 
structed the nonintersecting Brownian motions for any 
finite number A^ of particles all starting from the origin. 
There the Schur function plays an important role to rep- 
resent the transition probability density. We have shown 
that the spatial distribution of particles depends not only 
the observation time t but also on the time interval T in 



Assume 1 < n < m and let z(j, k), 1 < j < n,l < k < 
m, be indeterminates and Z(z(j, k)) be the nx m matrix 
with (j, A:) -element z(j,k). We consider the sum of all 
minors of Z with a given size r. That is, we define 

d(ai,- ■ ■ ,a r ) = V" det (z(aj,b k )) 

* — * Kj.k<r 
l<b 1 <b 2 <-<b r <m 

for r = 1, 2, • • • , m. Okada proved the following equalities 
known as the minor- summation formula [|35[ ; if r is odd, 
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then 



d(oi, ■ •• ,a r ) 

t 

-d(oi) 



Pf 



d(ai) d(a 2 ) 
d(ai,a 2 ) 



-d(a 2 ) -d(ai,a 2 ) 
\— d(a r ) —d(ai,a r ) —d(a 2 ,a r ) 



d(a r ) \ 
d(ai, a r N 
d(a 2 , a r 

o y 



if r is even, then 

d(ai, ■■■ ,a r ) = 

/ d(ai,a 2 ) d(oi,o 3 ) 

-d(a\,a 2 ) d(a 2 ,a 3 ) 

Pf -d(ai,a 3 ) -d(a 2 ,a 3 ) 

\— d(ai,a r ) —d(a 2 ,a r ) ~d(a 3 ,a r ) 



Let 

d(xj-) = 

d(xj,x k ) = 
and set 



^ «(xj,y(6i)), 

l<bi<m(M,5) 

y- z(xj,y(bx)) z(xj,y(b 2 )) 

i**^r? z ( x ^v( b i)) z(x k ,y{b 2 )) 

l<bi<b2<'m{M,5) 



d(ai,a r )\ 
d(a 2 , a r N 
d(a 3 , a r 

J 



d(xj, Xk) if 1 < j < k < n, 

—d(xj, Xk) if 1 < k < j < n, 

Sjk({%t}) = \ d(xj) if 1 < j < n, k = n + I, 

—d{xj) if j = n + 1, 1 < k < n, 

if 1 <j = fc < n + 1. 



t\ t , , . , , ■ • r rvi j i r Then Okada's formula gives 

Mow we give the integration version of Okada s for- to 

mula. Let z(x,y) be a square integrable continuous func- 
tion of real variables x, y. Then 



/ 



d n y det (z(Xj,y k )) 

oo<y 1 <-<y n <oo l<J,fc<n 



d n y det (z(xj,y k )) 

oo<y 1 <-<y n <oo l<],k<n 
Ph<j<k<n(F jk ({ Xe })), (A.l) 



lim lim S n Ph< j<k <n(S jk {{x e })). (A.2) 

M— >oo 5— *0 



where 
and 

^jfc({^}) 
with 

^js \Xj , X/c) 



n 

n + 1 



\P^j •> ^k) 
Iz \ i %k ) 

Iz(x 3 ) 
-Iz(x k ) 





z{x 3 ,y)dy, 



-oa<y 1 <y 2 <oc 



if n is even 
if n is odd, 



if 1 < j < k < n, 

if 1 < k < j < n, 

h 1 < j < n, k = n + 1, 

if j = n + 1, 1 < k < n, 

if l<j=fc<n+l 



z{x 3 ,yi) z{xj,y 2 ) 
z(x k ,yi) z(x k ,y 2 ) 



Since the Pfafhan in (A.2) is a finite summation of 
n/2 products of d(xj, x k Ys if n is even, and it is a finite 
summation of (n— l)/2 products of d(xj, x k ys multiplied 
by d(xt) if n is odd, we may have 



A, det. {z{xj,y k )) 

oo<y 1 <-<y n <oo l<J.fe<n 



Pfi<i<fe<ft lim lim <5««.*) S ifc ({a:/}) , 



where 

a(j, k) — 2 for 1 < j < fc < n, 
a(j, n+1) = 1 for 1 < j < n. 



dyidy 2 Since z ( a; ' 2/) ^ s assume d to be square integrable and con- 
tinuous. 



The p roof is the following. We write the integral in 
LHS of (A.l) as a limit of summation; 



d"y det {z{x 3 ,y k )) 

-oo<j/i<-<i/„<oc i S3-i k Sn 



= lim 



d n y det (z(xj,y k )) 



M-»oo J-M/2< yi <---<y n <M/2 l<j,k<n 



lim lim (55^+1 ({x £ }) = isfo), 

M— >oo 5— >0 

lim lim 5 2 5 ife ({^}) = I z (xj,x k ), 



for 1 < j,k < n. Then the proof is completed. 
By elementary calculation we can show that 



lim lim S n ^ det (z(xj,y(b k ))), 

M—>oc d— *0 * — * Ki,/c<n 
l<6l<---<b N <m(M,5) ~ ~ 



where m(M,S) = [M/6], the greatest integer not greater 
than M/5, and 



h{x) = 1, 
h{x,y) = 



=Erf 



2^ 



M -5 



i(M, 5) -I { m(M, S) - 1 



M-S 



M 



for z(x, y) 



/2?Tt 



-(x-y) 2 /2t 



Applying (A.l) with n 



N the expression (ETj) is obtained from (t 
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